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Abstract 
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1 Introduction 

Much attention has recently been devoted to the perturbative "hard", or BFKL, pomeron 
[1], especially in relation to the study of the small x behaviour of the deep inelastic scattering 
structure functions (see a recent review in [2] ) . In application to soft phenomena, the value of 
the pomeron intercept is of principal importance. For the BFKL pomeron it is considerably 
above unity: 

o:bfkl(0) = 1 - (3a s /7r)£b, 

where the "energy" Eq is equal —4 In 2, and a s is the (fixed) QCD coupling constant [1]. 
However the BFKL pomeron is only the simplest of the family of pomerons (with positive sig- 
nature) and odderons (with negative signature) formed in a system of n interacting reggeized 
gluons. Should some of them also result supercritical (with the intercept above unity), the 
study of the high energy behaviour of the QCD cross-sections would require summation of 
an arbitrary number of exchanges of all these supercritical pomerons and odderons. 

It is extraordinarily difficult to obtain an explicit solution or an exact energy value for 
n > 2 interacting gluons (see, however, some ideas in [3]). For n = 3 (odderon) application 
of conformal symmetry allows to reduce calculations to an one dimensional problem [4]. 
Variational treatment with a simple trial wave function then gives an intercept also above 
unity, although lower than for the pomeron [5]: 

a o dd(0) = 1 - (3a s /ir)E odd , E odd < -0.37. 

For larger number n of pomerons different crude approximations lead to wildly different 
estimates for the intercept. Developing the interaction in powers of the conformal group 
Casimir operator an intercept linearly rising with n is obtained [6] . On the other hand, the 
Hartree-Fock approximation, presumably valid for large n, gives intercepts less than one and 
linearly falling with n [7]. In this approximation the energy is e rising linearly with n: 

E n = (n/2) * 0.959... 

Evidently negative energies, as witnessed by the pomeron and odderon, result from cor- 
relations. In absence of explicit solutions for n > 2 gluons and in view of great technical 
difficulties in working with conformally invariant wave functions, evident already for the 
odderon case, the only realistic approach for n > 3 gluons seems to be a direct variational 
one, based on some complete and simple basis of functions. Taking a finite number N of 
these, one then computes a finite energy matrix. The intercept (minus one) is found as its 
smallest eigenvalue (ground state) with an opposite sign. With N growing, the ground state 
energy goes down, so that one always obtains its upper limit. As N — ► oo one is sure to 
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obtain an exact value, provided the basic functions form a complete set. One can have some 
idea of the precision for a given N by calculating the energy for the two-gluon case where its 
exact value is known. 

In this paper we give an outline of such a variational method aimed at calculating of the 
intercept of the system of an arbitrary number of reggeized gluons in a symmetric configu- 
ration both in the ordinary space and in the colour space. The basic functions were chosen 
as harmonic oscillator functions in the variable z = lnr 2 (radials) multiplied by azymuthal 
functions for given angular momenta / = 0, ±1,±2, .... They retain some of the conformal 
invariance corresponding to the substitution r — > 1/r. To see the convergence, we studied 
the two-gluon (the BFKL pomeron) and three-gluon (odderon) cases, where the exact en- 
ergy value for the former and its upper limit for the latter are known. It turns out that the 
convergence is rather slow. In the present series of calculations up to 201 basic states have 
been included for the pomeron and up to 1335 basic states for the odderon. With these basic 
states the resulting upper limit achieved for the two-gluon ground state energy (the BFKL 
pomeron) is —1.032, which means ~ 53% of the exact correlation energy. The upper limit 
obtained for the odderon ground state energy is +0.331. This limit lies considerably above 
the one obtained in [5] with a conformally invariant ansatz. From the known BFKL intercept 
we can study the dependence of the calculated energy on the number of basic states taken 
into account. Applying a similar fit to the odderon case and extrapolating to infinite number 
of states included gives an estimate 

Eodd < -0.3 - -0.6, 

in agreement with [5]. 

2 Variational calculation of the ground state energy for n 
reggeized gluons 

As shown in [8] the transverse space and colour wave function ip of n reggeized gluons in 
a colourless state satisfies a Schrodinger-like equation 

n 

H^ = El[ P ^, (1) 

1=1 

where pi is the momentum of the i-ih gluon. The Hamiltonian H is given by a sum of pair 
terms 

H = -{\/Q)Y,TiT k H ik . (2) 

i<k 
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Here Tj is the colour vector of the i-th gluon. In a colourless state 

n 

i=i 

The pair Hamiltonian acts on the wave function according to 

n n 

H ik 4> = H P 2 j (\n P y k + 4C)^+ J] V ){p\ ln(ri/4)pg + (i <_> fc))^ + 2(p i +p fc )V(ri fc = 0), 

(4) 

where rjfc = n — is the (transversal) distance between the gluons; C is the Euler constant. 
The energy eigenvalue E is related to the intercept by 

a(0) = 1 - (3a s /vr)£, (5) 

so that the rightmost singularity in the complex angular momentum plane corresponds to 
the ground state. As mentioned, for n = 2 the solution of (1) is the BFKL pomeron with 

E = -4 In 2, a BFKL {Q) = 1 + (a s /vr)12 In 2. (6) 

No explicit solution has been found for n > 2. 

The solution of (1) may evidently be found by a variational approach, searching the 
minimum value of the functional 

*= ( Hd 2 Pi rH^ = j2^k, (7) 

J i<k 

with the normalization condition 

/n 
n d w*n^=i- (8) 
i=i 

The BFKL pomeron and odderon states are symmetric both in colour and space variables. 
For more gluons solutions do not presumably possess such simple symmetry properties, since 
Eq. (1) mixes colour and space variables. Among various solutions there certainly are those 
which correspond to a number of BFKL pomerons weakly interacting with each other (mul- 
tipomeron cuts in the old Regge-pomeron theory). These solutions are basically symmetric 
neither in space nor in colour variables. We are interested in solutions of a different type, 
which do not reduce to BFKL pomerons and represent new pomerons with a higher intrinsic 
colour. While exact symmetry properties of these new states are not evident, one expects 
that in the variational approach they can well be approximated by states symmetric both in 
colour and space variables, similar to the simplest BFKL pomeron and odderon states. In 
any case, the energy value obtained with such a choice will give an upper limit for the exact 
one. According to the colour cluster separation property discussed in [7], these new pomerons 
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(and odderons) will also appear as subsystems for still larger number of gluons, that is, will 
take part in multipomeron exchanges. Such a picture was phenomenologically introduced in 
a model of colour string fusion [9], where it was supposed that supercritical pomerons exist 
for arbitrary high intrinsic colour. 

As noticed in [4], for a wave function symmetric both in colour and space variables, 
instead of the full functional <E> one can use any of its pair terms ^ in (7) with i and k fixed. 
Moreover, for a symmetric colour wave function, 

< T,{T k >= (2/n(n - 1)) ]T T{T k = -3/(n - 1), (9) 

i<k 

where (3) and T 2 = 3 have been used. As a result, the energy E may be found via the 
minimal value of the functional in only spatial variables 

£ = (1/2) J f[ d 2 Pi rH 12 i>, (10) 
i=i 

where the Hamiltonian H\2 is defined by (4) with ik = 12 and the function ip should satisfy 
(8) . The energy of the whole system of n gluons is determined by the minimal value e of £ 
according to 

E n = (n/2)e n . (11) 

Actually the operator in £ does not depend on n. The dependence on n enters only from 
extra arguments in -ip through the requirement of the symmetry in all arguments. The space 
of trial wave functions then gets smaller with rising n, wherefrom one obtains 

e«+i > e n (12) 

and in particular e n > e 2 - This gives a crude lower limit on the energy and an upper one on 
the intercept for n gluons in a symmetric state 

E n > (n/2)E . (13) 

This limit was obtained in [4] for the odderon n = 3. 

For a large number of reggeized gluons one can hopefully apply the Hartree-Fock approx- 
imation and seek for the minimum of the functional £ on functions ip which factorize into 
a product of individual gluon wave functions. For the ground state, with tjj symmetric in 
all gluons, all individual gluon wave functions should evidently be the same. As calculated 
in [7], the obtained minimal Hartree-Fock value is e = 0.959..., so that the Hartree-Fock 
energy of the symmetric n gluon state is positive and grows linearly with n, as stated in the 
Introduction. It means that there is little hope to expect supercritical pomerons composed of 
a large number of gluons. It also means that supercritical pomerons are formed exclusively 
due to correlations in the wave function. 
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For a finite number n of gluons one can expand the symmetric spatial wave function in a 
sum of products of individual gluon functions: 



n 



ip(r 1 ,...,r n )= Coi,-,«n Tl^aM), (14) 

ai,...,a n i=l 

where the one- gluon functions ip a (ri) form a discrete complete set and are orthonormalized 
according to (8): 

J d 2 r7plp 2 tp a > = 5 a:a >. (15) 

The coefficients c aiv .. jCtn have to be symmetric in all a's by the requirement of the Bose 
symmetry and normalized according to 

E |G«i,-,aJ 2 = 1 (16) 

The two-gluon Hamiltonian #12 acts nontrivially only on the wave functions for the gluons 
number one and two. So the energy functional becomes 

£ ~ C oi 1 ,a2,a2,,...,a n c a' 1 ,a' 2 ,a 3 ,...,a n ^a 1 ,a2,a' 1 ,a' 2 -i (17) 



ai,a2,a' 1 ,a 2 ,a3,...,a„ 



where the matrix £ ai)a2jQ ' ja ' is the two-gluon energy in the basis formed by functions tfj^ . 
With this matrix known, the problem of minimization of the functional £ reduces to finding 
the minimal value of a cuadratic form, that is, the minimal eigenvalue of the matrix £ aija2jQ ^ )a ^ 
considered as a matrix in independent initial and final n-gluon states. The latter means 
that this matrix should be multiplied by unity matrices for the rest of the gluons and then 
symmetrized in all initial and final gluons. The procedure is quite straightforward, once 
the basic functions ip a are chosen. It however involves a numerical evaluation of the energy 
matrix elements and a diagonalization of the matrix, whose dimension is rapidly growing 
with the number of gluons and the basic functions taken into account. 



3 Two-gluon energy matrix for given angular momenta 

The first task in the calculation of the energy matrix is to separate radial and angular 
dependence. The basic functions depend on the azymuthal angle <p trivially: 

tp a (r) = tp k ,l( r ) ex P Ufa ( 18 ) 

where I = 0, ±1, ±2, ... and k = 0, 1, 2, ... enumerates the radial functions. Thus a = {k, 1} is 
in fact a pair of indices. In the following, instead of r, we shall use the variable z = lnr 2 in 
most cases. In terms of z and 4> 

= -(4/r 2 )(d 2 z + (1/4)^), (19) 
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so that acting on the function tpk,i{r) 

pV*,i(r) = -(4/r 2 )(d 2 - (l/4)Z 2 )^ M (r). (20) 

Wave functions with different values of the angular momentum are automatically orthogonal. 
For coinciding I the normalization condition for the radial functions takes the form 

J dzr k ,i(z)(-d 2 + (1/4)Z 2 )W,^) = (1/4tt)<W. (21) 
It reduces to the standard form for functions 

e*,j(*) = (0+IWiw(*)> (22) 

which evidently satisfy 

/ dzC Kl {z)^ :l {z) = (1/4tt)^. (23) 

In the following we assume that the radial functions are chosen to be real. 

With the angular dependence of the wave function explicitly given by (18), one can do 
the azymuthal integrals in the potential energy in a straightforward manner. Let aj = {ki, k} 
and take the transition between two gluon states ai,(X2 — ► «3, «4- Evidently the total angular 
momentum is conserved so that the energy matrix elements are zero unless h + h = h + 
According to (4) the potential energy consists of two parts, the first part U given by an 
essentially Coulomb interaction between the gluons and the second one Q given by a contact 
interaction, proportional to their total momentum squared. Let us begin with the Coulomb 
part U. Its two terms evidently give the same contribution due to the symmetry under the 
interchange of gluons 1 and 2. So we can take only one of them and drop the factor 1/2. 
Denote 

= (d 2 - (l/4)Z 2 )Y> M (z). (24) 

Then after doing the azymuthal integration we obtain the following matrix elements for the 
potential energy U 

= 16vr 2 J dz 1 dz2Va l (zi)ijJa 3 {zi)tpa 2 {z2)riaA z 2)Ui{z 1 ,Z2), (25) 

where I = \l± — Z3 1 = \fo — h\ is the angular momentum transfer and the function Ui(z±, z-i) is 
given by 

l/, = -(l/0exp(-(i/2)(|zi-Z2|), 1^0, (26) 

and 

Uq = maxjzi, z 2 }. (27) 
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The contact part Q involves gluonic wave functions taken at the same point, that is, with 
the same r and <f>. After performing the azymuthal integration and integrating once by parts 
in the variable z we obtain 

Q 'qi ,Q2;03 ,0-4, 

= 87T 2 I dz({d + (l/2)|Ji + hWa^aMd + (V2)|*3 + WasW- ( 2 §) 

One can somewhat simplify this expression by noting that 

(8 + (1/2) |Zi + hD^mlpai = £ai^a 2 + ^ai£a 2 + A 12 1pai^a 2 , (29) 

where 2A12 = |/i + h\ — \h \ — \h\ and similarly for the second factor in (28). Then finally 
= 8tt 2 J dz^ ai if> a2 +if> ai ^ a2 + A 12 ip ai ^2)^3^4+^3 ^4 + ^34^3^4)- (30) 

Passing to the kinetic energy given by the first terms in (4) we also note that the two 
terms give the same contribution so that we can take one of them (say, for the gluon 1) and 
drop the factor 1/2. For the transition of the first gluon a.\ — > 03 the conservation of angular 
momentum requires that l\ = I3. The kinetic energy is easily calculated in the momentum 
space. So we transform the basic functions to the momentum space according to 

Mp) = J(d 2 r/2Tr)^ a (r) exp(-ipr). (31) 

The azymuthal integration leads to 

^a(p) = (-i) l expilcf) f rdri> k j(z)Ji(pr), (32) 

where Ji is the Bessel function. To do the integral over r it is convenient to introduce a 
Fourier transform of the function ip with respect to the variable z: 



4>k,l{z) = J {dv / \/ 2ir)4> kt i(v) expiuz. (33) 
Putting this representation in (33) and doing the r-integration we obtain 

W(P) = (2/p 2 ) exp il<f> J (dv/V2^)f k j(v)p- 2l », (34) 

with 

fkM = (-ifh 2w (\l\/2 + iu)^> Kl {u)Y{\l\/2 + iu)/Y{\l\/2 - iv). (35) 

With the gluon wave functions in the momentum space given by (34), both radial and 
azymuthal integration in p are easily done. The final matrix element of the kinetic energy T 
results as 

T ai , as = -4™ J duf* kuh (u)(d/du)f k3h (u) (36) 
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(recall that l\ = Is). The differentiation gives 

(d/dv)f k3:h (v) = / fc 3^ 1 ( l /)(2zln2 + 2iR e ^(|/i|/2 + ^) + (a/^)ln((|/|/2 + ^)^ M (^)))- (37) 
Correspondingly the kinetic energy matrix element separates into terms 

T&as =^ J ^/fe 1 AWA3,/i(^)(ln2 + Re^|/i|/2 + ^)) (38) 

and 

Tg>,« 3 = J ^((|/i|/2 + ^)^ 1 , il (^))*(a/^)((|/ 1 |/2 + ^)^, il (^)). (39) 

One notes that the function ((|/|/2 + iv)(t>k,i{v)) is nothing but the Fourier transform of 
Z,k,l{z) w hh respect to z. Correspondingly we denote it as 

(|Z|/2 + «/)<^) = &>)- (40) 
The part can then be written as 

Ti% 3 = -Am J du^ uh (uy(d/du)^ h (u)). (41) 

The orthonormalization property (23) transforms into the analogous property in the v 
space 

J dvt,l l (v)(i k ,, l (u) = {l/AK)5 kk ,. (42) 

Noting that f kl (v)f k ',i{v) = £| i( v )£,k' we observe that the term ln2 in (38) will add a 
constant 2 In 2 to the energy. Separating another constant term 2-0(1) we finally present the 
part in the final form 

T« Q3 = 2(ln2 + 0(l))5 Ql , Q3 +8^ J dvC ai {v)U{u){Re ^(\h\/2 + iv) - ^(1)). (43) 

The first constant term cancels an identical one in the initial Hamiltonian (4), so that we 
may forget about these constants and concentrate only on the resting nontrivial contributions. 
Using the representation 

roo 

ip{x)-ip{l)= dt(exp(-t) -exp(-xi))/(l - exp(-t)) (44) 
jo 

and the othornormalization property of the set £ Q we may cast in the form 

roo 

Ti2«s = 2 / W(expt - l))(6 ai a 3 ~ ex P (i(l - |Zi|/2))<fa ia3 (t)), (45) 

J 

where 

5aia 3 (*)= 47r / dv^{v)^, a3 {v)cosvt. (46) 
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Note that (44) is not valid for Rex = 0. Therefore this formula cannot be applied when the 
gluon orbital momentum is zero. In this case one may use 

tp(iv) + ip(—iv) = ip(l + iv) + ip(l — iu), 

which formally corresponds to changing the angular momentum to be equal to 2. 

As to the second part of the kinetic energy , we shall find out presently that it will be 
cancelled by a similar contribution coming from the monopole part of the Coulomb interaction 
for the angular momentum transfer equal to zero. 



4 Monopole part of the Coulomb interaction 

Most of the contributions to the energy presented in the previous section can hardly be 
further simplified and were used in the numerical calculations as they stand. The exception 
is the monopole part of the Coulomb interaction corresponding to (25) with I = (Eq. 
(27)). This part contains contributions which cancel the term in the kinetic energy 
and partially the contact interaction contribution for / = 0. The cancellation between the 
monopole Coulomb interaction and the kinetic term is responsible for the scale invariance 
of the energy. 

Explicitly the monopole term contribution is given by 

/oo rzi 
dz 1 ri ai (z 1 )ip ol3 (z 1 )z 1 / ^ a2 {z 2 )r] a4 {z 2 ) + (ai <-> a 4 ,a 2 <-> a 3 ). 
-OO J —CO 

(47) 

Here and in the following it is assumed that I = 0, that is, l\ = 1% and l 2 = h- Introduce a 
function 

Xa 2 ,a 4 (z)=[ dz'ip a2 (z')f] a4 (z'). (48) 

Once integrating by parts we find 

Xa 2 ,a 4 (z) = ^ a2 (z)£ a4 {z)- 6* 2 ,a 4 (z) , (49) 

where the function ^ a2i a A {z) with two indices, symmetric in these, is defined as 

ta 2 ,a 4 {z) = ( dz'£ a2 (z')£ a4 (z'). (50) 
J — oo 

As z -> oo, according to (23), £, a2y a 4 {z) -> (l/47r)<5 a2iQ4 , so that 

Xa 2 ,a 4 (co) = -(l/47T)(5 Q2iCe4 . 

Having this in mind we can rewrite (47) in the form 

/'OO 
dz(xa 3 , ai (z) + (l/47r)5 Ql 
,az) ZXa 2 ,a 4 
(z) + (c*i <-> 04, a 2 <-> 03)- (51) 
-00 
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Integrating by parts, the integral transforms into 

/oo 
dz{Xaa,ai(z) + (l/^)5 aija3 )(zx a2 ,a 4 (z) + Xa 2 ,a 4 )- (52) 
-oo 

The term coming from the product Xa3,ai^Xa2,a 4 cancels the contribution (cci <-> «4, «2 CK3) 
in (51) so that the monopole contribution becomes 

/oo 
03,01 (^)Xct 2 ,C*4 
(z) + (l/47r)<5 aiiQ3 
(*)))■ (53) 
-00 

Now the idea is to substitute the functions x in (53) by the symmetric functions £ using 
relation (49). Take the the first term in the integrand of (53). With (49) we obtain for it 

Xa3,aiXa2,«4 = ^0:3 £ai V'a 2 £a4 — V^Cai^o^,^ — £a3,ai^a2£a 4 + £0:3,0:1 £a 2 ,0-4 • 

Having in mind the subsequent symmetrization with respect to the interchange of gluons 1 
and 2, we can change a,\ <-> ai and 03 <-> 04 in the second term. Summed with the third 
term it then gives 

- £03,01(^6*4 + £02^04)- (54) 

Recall now that £ Q2 = (3 + (1/2) \h\)^a 2 an d similarly for £ a4 . Integration by parts allows to 
substitute (54) by 

(£o 3 £oi - |k| £03,01)^02^04- (55) 

So finally the first term in (53) leads to the following three contributions to the monopole 
Coulomb energy: 

f°° 

U^ :a2 . a3 ,a 4 = -!6vr 2 / cte£ Q3 , ai £ Q2iQ4 , (56) 
J —00 

/oo 
dz£ az , ai il> a2 il> ai (57) 
-00 

and 

roo 

^oi,02;o3,04 = ~ 167r / <MY> a3 £ ai Y> a2 £o 4 + £036*1^02 ^o 4 ) • (58) 
J — oo 

Of these terms the first is divergent in its present form. It will receive its meaning after 
adding a new contributions coming from the rest of the terms in (53). For that reason we 
have denoted it with a tilda. 

Now for the rest of the terms in (53). Changing the function x by £ according to (49) we 
have under the integral 

£o 2 ,o 4 + z £o 2 ,o 4 = ^026*4 ~ £o 2 ,o 4 + ^a 2 ( d ~ (V 2 ) \h |)£a 4 • 

Integration by parts transforms it into 

£ct 2 ,o 4 ^£ct 2 £a4' (59) 
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The first term can be combined with (56) to give the final part U^: 

roc 

tffiUas.cu = 167r2 / <feU,«4((l/47r)<S Ql)a3 - U.aJ- (60) 

J— oo 

Now the integral converges due to the property (23). Putting here the explicit form of the 
functions i ai ,a k and integrating over z we obtain the term L^ 1 ) in its definitive form: 

U ai,a 2 ;a 3 ,a 4 = j dz X dz 2 {z X ~ Z 2 )9(z 1 - Z 2 )£ ai (2l)fa s («l)Ca 2 (^)Ca 4 W (61) 

The second term in (59) gives the last contribution to the monopole energy 

/oo 
dzz^ a , 2 ^ a4 . (62) 
-oo 

This term cancels with the contribution to the kinetic energy. Indeed after the Fourier 
transformation to the v space according to (33), the factor z goes into idjdv. One can then see 
that (62) gives exactly the contribution T^ 2 \ Eq. (41), with an opposite sign and with gluons 
1 and 2 interchanged, which is of no importance because of the subsequent symmetrization. 

The term cancels with the part of the contact interaction Q, Eq. (30), which does 
not contain factors A: 

2 l°° 

Qai,a2\az,a4 ^ / dz(l(j a3 £ ai 'llj a2 £ ia4 + ^az^a-i^P ct2^ a-A + V ; aaV'aiCa2?a4 £«3^ 'aiCa^V '0:4) " 
J -co 

(63) 

Summed with this part gives 

/oo 
dz(4> ai €a 2 - CaiY'aaXCasYW + ^a 3 £a 4 )- (64) 
-00 

This expression is antisymmetric under the interchange of the gluons 1 and 2 and does not 
give any contribution to the energy. 

So finally the only contributions which remain in the interaction for zero angular momen- 
tum transfer are and the part of the contact interaction (30) which contains 
factors A. 

5 Oscillator basic functions 

A natural orthonormal discrete basis for z = lnr 2 varying from — oo to +oo is formed by 
the harmonic oscillator proper functions. Thus we choose functions £,k,i(z) independent of I 
and given by 

= c k H k (z)exp(-z 2 /2), (65) 

where H k are the Hermite polinomals and c k are determined by the normalization condition 
(23) to be 

c 2 k = l/(4vr 3 / 2 2 fe fc!). (66) 
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The Fourier transformation to the v space gives 

& W = (-i) k c k H k (u) exp(-i/ 2 /2). (67) 
In the coordinate space the function r] k j(z) is obtained from £ by differentiation: 

r)kjl (z) = (d- (l/2)|Z|)&(z). (68) 
Using the properties of iJ^ we get 

= 2fc(c fc /c fc _ 1 )£ fe _ 1 (z) - (z + (1/2)|*|)&(*). (69) 
The function i/j k ,l is obtained from as a solution of the differential equation 

U^) = {d + {l/2)\l\)^ l {z), (70) 
with a boundary condition ip k: i(—oo) = 0. It is given by an integral 

y, M (z)=f ^ / a(^)ex P (-|Z|(z-/)/2). (71) 

J— oo 

For = 0, 1 we find from (71): 

^o,l = c exp((u; 2 - z 2 )/2)(l - ${w/y/2)), (72) 

i; hl (z) = 2(z + w){ Cl /co)^o,l( z ) ~ 2c i exp(-z 2 /2), (73) 

where w = \l\/2 — z and <& = erf is the error function integral. For k > 1 the functions tp can 
be found by a recurrency relation that follows from (71) upon integrating by parts: 

i>k+i,l = \l\(c k +i/c k )ipk,l + 2fc(cfc + i/c fe _i)V'fc-i,i - 2(c fe+ i/c fc )^ fe . (74) 

With this set of functions the potential part of the energy was calculated numerically 
As to the kinetic part, the function g entering (45) can be found analytically. For transition 
k, I — > k',l it is equal to zero if k + k' is odd. For even k + k' = 2s and k > k' 

k' 

g kk ,(t) = 4vr 3 / 2 (-l)<W exp(-i 2 /4) £ 2Pp!C*C* (-t 2 )'-* , (75) 

p=0 

where 2d = k — k'. 

We finally note that the presented set of basic functions can be trivially generalized to 
include a scaling factor a in the z space by choosing as basic functions 

fW = ^M. (76) 

Other functions can then easily be found as 

£\v) = (l/y/a)t k {u/a), (77) 
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%i( z ) = aVaVk,i/a(az) 



(78) 



and finally 



^k\ z ) = 0-Va)ipk,l/a(az). 



(79) 



Putting these expressions into the obtained formulas for the energy and rescaling the variables 
z or v we find that the new value of the potential energy is obtained by dividing all angular 
momenta by a and dividing the resulting energy by a. As to the kinetic energy the change 
reduces to substituting at instead of t in the function g, Eq. (45). 

6 Results for n = 2,3 and conclusions 

We applied this formalism to the cases n = 2 (the BFKL pomeron) and n = 3 (the 
odderon) to study its convergence. 

The energy matrix £ aita2t0 / a ' has been calculated numerically for a chosen set of basic 
functions described in the previous section. After proper symmetrization in two- or three- 
gluon states its lowest eigenvalue has been determined, which gives an upper limit on the 
exact pomeron or odderon energy according to Eq. (11). To study the minimal energy only 
states with the total angular momentum equal to zero have been included. 

The selected set of basic one-gluon functions is characterized by the maximal value of the 
angular momentum included l max and numbers of radial functions included for each wave. As 
calculations show, best results are obtained when one raises l max and the number of radials 
in all waves simultaneously. So we present here the results for the case when the number 
of radials r is the same for all angular momenta and is equal to the number of angular 
momenta included r = l max + 1. Such a set of functions is thus characterized by a single 
parameter r. With a growth of r the number of states N rises very rapidly. For two gluons 
N 2 = r(r 2 - r/2 - 1/2) and so rises as r . For three gluons the rise is still steeper. 

In the present calculations the number r was limited by 6 for two gluons and by 5 for 
three gluons. Correspondingly the total number of basic states included was taken up to 201 
for two gluons and 1335 for three gluons. 

The results of the calculations of the ground state energies of the BFKL pomeron e 2 and 
the odderon €3 are presented in the Table for different values of r. One observes that the 
obtained energies are still rather far from the exact value for the pomeron and the upper 
limit for the odderon obtained in [5]. Thus the convergence of the method is rather slow. 
The Table also reveals that the odderon energy is essentially larger than the pomeron one for 
a given r. So our results confirm that, in all probability, the odderon intercept is lower than 
that of the BFKL pomeron. 
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To be more quantitative one can estimate the precision of the variational results by com- 
paring the calculated correlation energy with its exact value known for the BFKL pomeron. 
Subtracting Eq from the Hartree-Fock energy we find e!f r = 3.732. For the maximal value 
r = 6 the calculated correlation energy is 1.991, which constitutes ~ 53% of the exact value. 
With r = 5 the correlation energy is 1.871, that is ~ 50% of the exact value. If one boldly 
assumes that for the odderon the correlation energy calculated with r = 5 also constitutes 
50% of the exact value, then one finds this exact value to be 1.256 and the absolute energy 
€3 = —0.3, that is, E oc ^ = —0.45, which does not contradict [5]. 

More sophisticated estimates can be attempted by studying the dependence of the cal- 
culated energy on r and extrapolating for higher values of r. We have chosen a 4-parameter 
fit 

e(r) = e(oo) + aexp(— a In 13 r). (80) 

For the BFKL pomeron e(oo) is known. It turns out that the values of e^r) given in the 
Table are well described by (80) with the choice of parameters 

a = 6.066, a = 0.925, (3 = 0.515. 

According to this fit further improvement of the value for ei requires very high values of r. 
E.g., to achieve e2(r) = —2.0 one has to raise r to ~ 100. The corresponding numbers of 
two-gluon states are enormous and hardly possible to include. 

The analogous fit for the three-gluon case determined from energy values for r = 2, 3, 4, 5 
has the parameters 

e(oo) = -0.389, a = 2.472, a = 1.04, (3 = 0.36. 

However this set of parameters is rather unstable: a small change in the energy values used 
causes rather large changes in the value of e(oo). If we take for (3 = 1/2, as evidently favoured 
by the two gluon case, then we obtain 

e(oo) = -0.195, a = 1.772 a = 0.957, (3 = 0.5, 

with the value €3(2) = 0.603 slightly smaller than the calculated one 0.605 although within 
the calculational errors of the order of ±0.003. In view of this we can only give a crude 
estimation for the upper limit for the odderon energy from our calculations, 

e 3 < —0.2 -. — 0.4, 

which according to (11) translates into the estimate for E ^ 



E odd < -0.3 - -0.6, 
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in agreement with [5]. 

To conclude, the calculations for two and three gluons show that the developed method can 
be applied to investigate the intercept of symmetric configurations although the convergence 
is slow, evidently, due to a very singular character of the gluonic wave function. It seems 
realistic to obtain around 50% of the correlation energy with this approach, which may serve 
to estimate intercepts for multigluon configurations relative to the BFKL one. We are trying 
to achieve better results by including more basic functions and also selecting configurations 
which give the dominant contribution. Calculations for n = 4 gluons are also in progress. 
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Table 

Calculated values of the ground state energy per gluon (multiplied by 2, Eq. 
(11)) for the pomeron (62) and odderon (€3) with different numbers r of radial 
functions and angular momenta included. 



r 


£2 


£3 


1 


0.968 


0.968 


2 


0.022 


0.605 


3 


-0.475 


0.454 


4 


-0.743 


0.379 


5 


-0.912 


0.331 


6 


-1.032 










00 


-2.773 





